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Abstract The fundamental description of both structural properties and reactions 
of light nuclei in terms of constituent protons and neutrons interacting through 
nucleon-nucleon and three-nucleon forces is a long-sought goal of nuclear theory. 
I will briefly present a promising technique, built upon the ab initio no-core shell 
model, which emerged recently as a candidate to reach such a goal: the no-core shell 
model/resonating-group method. This approach, capable of describing simultane- 
ously both bound and scattering states in light nuclei, complements a microscopic 
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cluster technique with the use of two-nucleon realistic interactions, and a micro- 
scopic and consistent description of the nucleon clusters. I will discuss applications 
to light nuclei binary scattering processes and fusion reactions that power stars 
and Earth based fusion facilities, such as the deuterium-^He fusion, and outline 
the progress toward the inclusion of the three-nucleon force into the formalism and 
the treatment of three-body clusters. 

PACS 24.10.-i • 24.10.Cn • 21.60.De • 25.40.Lw • 25.45.-z • 21.45.Ff 



1 Introduction 

Nuclei are aggregates of protons and neutrons interacting through forces arising 
from the underlying theory of quantum chromodynamics. Understanding how the 
strong force binds nucleons together in nuclei is fundamental to explain the very 
existence of the universe. Indeed, the mutual interactions between nucleons led 
to the formation of the lightest nuclei a few minutes after the Big Bang, and the 
following nuclear processes, producing heavier nuclei during stellar evolution and 
in violent events like supernovae, have been crucial in shaping the world we leave 
in. Therefore, one of the central goals of nuclear physics is to come to a basic 
understanding of the structure and dynamics of nuclei. The ab initio (i.e. from 
first principles) no-core shell model/resonating group method (NCSM/RGM) [ll 
[2] is a theoretical technique that attempts to achieve such a goal for light nuclei. 



2 Ab initio NCSM/RGM 

In the ab initio NCSM/RGM approach the many-body wave function. 



is expanded over a set of translational-invariant cluster basis states describing two 
nuclei (a target and a projectile composed of ^4— a and a < A nucleons, respectively) 
whose centers of mass are separated by the relative coordinate VA^^.a s-nd that are 
traveling in a '^^Ij wave or relative motion (with s the channel spin, I the relative 
momentum, and J the total angular momentum of the system): 

l^ur ) = \{\A-aaiI-i^^Ti)\aa2l2T2)) Yi[rA-a,a)\ ^• 

L J TA-a.a 

(2) 

Here, the antisymmetric wave functions lA—aail^'-Ti) and |aa2^2^^^2) are eigen- 
states of the {A — a)- and a-nucleon intrinsic Hamiltonians, respectively, as ob- 
tained within the NCSM approach [3] and are characterized by the spin-parity, 
isospin and energy labels I^\Ti, and o^, respectively, where i = 1,2. Additional 
quantum numbers labeling these RGM-inspired continuous basis states are parity 
TT = 7ri7r2(— 1)^ and total isospin T. In our notation, all these quantum numbers 
are grouped into a cumulative index v = {A—aail^^Ti; aa2l2^T2; s£}. The Pauli 
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principle is enforced by introducing the appropriate inter-cluster antisymmetrizer, 
schematically 

i. = y(^ (i+j:i-rp], (3) 

where the sum runs over all possible permutations of nucleons P different from 
the identical one that can be carried out between the two different clusters, and p 
is the number of interchanges characterizing them. Finally, the continuous linear 
variational amplitudes xt if) are determined by solving the orthogonalized RGM 
equations: 

E / drV-[Ar-^HAA-^].V(^,0^^ =i^^^, (4) 

v' 

where M^^T" {r^r'^ and T-L^^T {r^r'), commonly referred to as integration kernels, 
are respectively the overlap (or norm) and Hamiltonian matrix elements over the 
antisymmetrized basis i.e.: 

^u'7{r',r) = ($t:j\A,,A. \'l>i7) , ni:j{r',r) = (<l>i:j\ A,,HA, \<l>rj) . 

(5) 

Here, H is the microscopic ^4— nucleon Hamiltonian and E is the total energy 
in the center of mass (cm.) frame. For a detailed explanation of how norm and 
Hamiltonian kernels are obtained from the underlying nuclear interaction and 
the NCSM eigenvectors of target and projectile we refer the interested reader to 
Refs. g] and 0]. 

3 Applications 

Applications of the NCSM/RGM approach to the description of nucleon- and 
deuteron-nucleus type of collisions based on two-nucleon (NN) realistic interac- 
tions have already led to very promising results [1, 2, 5jl31[6l[7j. In most instances, 
we employed similarity-renormalization-group (SRG) [8l[9] evolved potentials, and 
in particular, those obtained from the chiral N'^LO [TD] NN interaction. Here we 
briefly review two of such applications, the calculation of the ''Be(p, 7)®B radiative 
capture [6], and the study of the ^H(d, n)*He and '^He(d,p)^He fusion reactions [7]. 

3.1 The ^Be(p,7)^B radiative capture 

Recently, we have performed the first ab initio many-body calculation of the 
^Be(p, 7)*B radiative capture [6], the final step in the nucleosynthetic chain leading 
to ^B and one of the main inputs of the Standard Solar Model. This calculation 
was carried out in a model space spanned by p-^Be channel states including the 
five lowest eigenstates of ^Be (the I ground and the ^ ,| , and first and second 

I excited states) in an Nra^^ = 10 NCSM basis, and employed the SRG-N^LO 
A''A'^ interaction with A = 1.86 fm^^, where A denotes the SRG evolution param- 
eter [8] . We first solved Eq. Q with bound-state boundary conditions to find the 
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Fig. 1 Calculated ''Be(p, 7)*B S-factor as a 
function of the energy in the center of mass 
compared to data. Only El transition were 
considered in the calculation. 
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Fig. 2 Convergence of the ^Be(p, 7)*B S- 
factor as a function of the number of 
^Be eigenstates included in the calculation 
(shown in the legend together with the cor- 
responding separation energy). 



bound state of ^B, and then with scattering boundary conditions to find the p-^Be 
scattering wave functions. Former and latter wave functions were later used to 
calculate the capture cross section, which, at solar energies, is dominated by non- 
resonant E\ transitions from p-^Be 5*- and D- waves into the weakly-bound ground 
state of *B. All stages of the calculation were based on the same HO frequency of 
TiQ = 18 MeV, which minimizes the g.s. energy of ^Be. 

At Nmax = 10, the largest model space achievable for the present calculation 
within the full NCSM basis, the '^Be g.s. energy is very close to convergence as 
indicated by a fairly flat frequency dependence in the range 16 < Tif2 < 20 MeV, 
and the vicinity to the A'^max = 12 result obtained within the importance-truncated 
NCSM [111112] . With the chosen value A = 1.86 fm^^ for the SRG evolution of the 
N^LO NN interaction, we obtain a single 2"*" bound state for ®B with a separation 
energy of 136 keV, which is quite close to the observed one (137 keV). This is 
very important for the description of the low-energy behavior of the ^Be(p, 7)^B 
astrophysical S-factor, known as S17. While for a complete ab initio calculation 
one should include also the three-nucleon (AAA) interaction induced by the SRG 
evolution of the AA potential as well as the SRG-evolved attractive initial chiral 
AAA force, we note that in the yl-range ~ 1.8-2.1 fm-\ and, in very light nuclei, 
the former interaction is repulsive and the two contributions cancel each other to 
a good extent p'SlfTi] . 

Figure 1 compares the resulting S17 astrophysical factor with several experi- 
mental data sets. Energy dependence and absolute magnitude follow closely the 
trend of the indirect Coulomb breakup measurements of Shiimann et al. [15| , while 
somewhat underestimating the direct data of Junghans et al. [16]. The resonance 
due to the All capture, particularly evident in these and Filippone's data and 
missing in our results, does not contribute to a theoretical calculation outside 
of the narrow ^B 1"*" resonance and is negligible at astrophysical energies [17] . 
The Ml operator, for which any dependence upon two-body currents needs to 
be included explicitly, poses more uncertainties than the Siegert's El operator. 
In addition, the treatment of this operator within the NCSM/RGM approach is 
slightly complicated by the contributions coming from the core (^Be) part of the 
wave function. Nevertheless, we plan to calculate its contribution in the future. 




Fig. 3 Calculated S-factor of the 
'^He((i, p)''He reaction compared to ex- 
perimental data. Convergence with the 
number of deuterium pseudostates in the 
^Si-^Di (d*) and 3Z)2 (d'*) channels. 



Fig. 4 Calculated 3H(d,n)4He S-factor 
compared to experimental data. Conver- 
gence with Afjnax obtained for the SRG- 
N^LO NN potential with A = 1.45 fm^^ 
at hS7 = 14 MeV. 



Our calculated 5*17(0) = 19.4(7) MeV b is on the lower side, but consistent with 
the latest evaluation 20.8 ± 0.7(expt)±1.4(theory) ^7\. The 0.7 eV b uncertainty 
was estimated by studying the dependence of the S-factor on the harmonic oscil- 
lator (HO) basis size A^'max as well as the influence of higher-energy excited states 
of the ^Be target. More precisely, we performed calculations up to A^max = 12 
within the importance-truncation NCSM scheme [111112] including (due to com- 
putational limitations) only the first three eigenstates of ''Be. The Nmax = 10 
and 12 S-factors are very close. In addition, the convergence in the number of ^Be 
states was explored by means of calculations including up to 8 ^Be eigenstates in a 
A'^max = 8 basis (larger A'^max values are currently out of reach with more then five 
^Be states) . This last set of calculations is presented in Fig. [2] from which it ap- 
pears that, apart from the two | states, the only other state to have a significant 
impact on the 5*17 is the second ^ , the inclusion of which affects the separation 
energy and contributes somewhat to the fiattening of the S-factor around 1.5 MeV. 
For these last set of calculations we used SRG-N'^LO interactions obtained with 
different A values with the intent to match closely the experimental separation 
energy in each of the largest model spaces. Based on this analysis, we conclude 
that the use of an A'max = 10 HO model space is justified and the restriction to 
five ^Be eigenstates is quite reasonable. 



3.2 The ^H(d, n)^He and ^He(d,p)^He fusion reactions 

In the following we present the first ab initio many-body calculations of '^H(d, n)'''He 
and '^He(d,p)^He fusion reactions [7] starting from the SRG-N'^LO A'' A' interac- 
tion with A = 1.5 fm~^, for which we reproduce the experimental Q-value of both 
reactions within 1%. These reactions have important implications first and fore- 
most for fusion energy generation, but also for nuclear astrophysics, and atomic 
physics. Indeed, the deuterium-tritium fusion is the easiest reaction to achieve on 
earth and is pursued by research facilities directed at reaching fusion power. Both 
'^H(d, n)'*He and ^He(d,p)'*He affect the predictions of Big Bang nucleosynthesis 
for light-nucleus abundances. In addition, the deuterium-^He fusion is also an ob- 
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ject of interest for atomic physics, due to the substantial electron-screening effects 
presented by this reaction. 

The model spaces adopted are characterized by HO model basis sizes up to 
A^max = 13 with a frequency of TiQ = 14 MeV and channel bases including n- 
^He (p-'*He), d-^H (d-^He), d*-^R {d*-^Be) and d'*-^H {d'*-^iie) binary cluster 
states. Here, d* and d'* denote ^Si-^Di and deuterium excited pseudostates, 
respectively, and the ^H (^He) and ^He nuclei are in their ground state. 

The results obtained for the '^He(d,p)^He S- factor are presented in Figure [3l 
The deuteron deformation and its virtual breakup, approximated by means of d 
pseudostates, play a crucial role in reproducing the observed magnitude of the 
S-factor. Convergence is reached for 9d* + 5d'* . The typical dependence upon the 
HO basis sizes adopted is illustrated by the '^H(d, n)'*He results of Fig. [H The 
convergence is satisfactory and we expect that an A^max = 15 calculation, which 
is currently out of reach, would not yield significantly different results. While the 
experimental position of the '^He{d,p)^He S-factor is reproduced within few tens of 
keV and we find an overall fair agreement with experiment (if we exclude the region 
at very low energy, where the accelerator data are enhanced by laboratory electron 
screening), the ^H(d, n)^He S-factor is not described as well with A = 1.5 fm~^. 
Due to its very low activation energy, the ■^H(d, n)^He S-factor, particularly the 
position and height of its peak, is extremely sensitive to higher-order effects in the 
nuclear interaction, such as the NNN force (not yet included in the calculation) 
and missing isospin-breaking effects in the integration kernels (which are obtained 
in the isospin formalism). With a very small change in the value of the SRG 
evolution parameter we can compensate for these missing higher-order effects in 
the interaction and reproduce the position of the '^H(d, n)^He S-factor. This led to 
the theoretical S-factor of Fig. [3] (obtained for A = 1.45 fm~^), that is in overall 
better agreement with data, although it presents a slightly narrower and somewhat 
overestimated peak. This calculation would suggest that some electron-screening 
enhancement could also be present in the ^H(d, n)*He measured S-factor below 
10 keV cm. energy. However, these results cannot be considered conclusive until 
more accurate calculations using a complete nuclear interaction (that includes the 
NNN force) are performed. Work in this direction is under way. 



4 Recent developments 

Here we outline some of our more recent efforts in the development of the NCSM/RGM 
approach, namely the progress toward the inclusion of the three-nucleon force into 
the formalism, and the treatment of three-body clusters and their dynamics. 



4.1 Scattering and three-nucleon force 

In past applications for light-ion reactions, we omitted the NNN interaction in- 
duced by the SRG renormalization of the A^'A'' potential as well as the initial chiral 
NNN force. By neglecting induced forces, we introduced a dependence on the SRG 
parameter A, which was then chosen so that the particle separation energies were 
well reproduced. For low-energy thermonuclear reactions, this is a dominant effect. 
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Fig. 5 Convergence with respect to the HO 
basis size Nmax = 11 at hO = 20 MeV of the 
n-|-*He{g.s.) phase shifts obtained for the 
SRG-(N3lO NN + N^LO NNN) interaction 
with A = 2.0 fm-l. The label A^ASmax = 14 
refers to the HO size of the NNN matrix el- 
ements. 




Ek,n [MeV] 

Fig. 6 Calculated n-|-*He(g.s.) phase shifts 
for SRG-N^LO NN-only (dots), NN-I-NNN- 
induced (dashed line) and SRG-(N3L0 NN 
+ N^LO NNN) interactions (solid hue) with 
A = 2.0 fm~^ obtained within the HO ba- 
sis size Afjnax = 11 and frequency hfi = 20 
MeV. See also the caption of Fig. [5] 



and overall such a procedure led to (never obtained before) realistic results. How- 
ever, a truly accurate ab initio description demands the inclusion of both induced 
and chiral NNN interactions. 

While the inclusion of the TViVA^ force into the NCSM/RGM formalism is 
conceptually straightforward, in practice it poses major challenges having to deal 
with: 1) the large number of NNN matrix elements, which makes it essential to 
work within the JT-coupled scheme; and 2) the appearance of kernels depending 
on the three-body densities of the target already for nucleon-nucleus processes, 
which demands new efhcient computing strategies for applications with p-shell 
targets to be possible. Figures [5] and |6] present initial results for *He(n, n)*He 
scattering phase shifts with inclusion of the NNN force. The use of SRG-evolved 
interactions facilitates the convergence of the calculation within A'^max ~ 11- At 
the same time. Fig. |6] highlights the influence of induced and initial components 
of the NNN force on the resonant phase shifts. The largest splitting between 
^-iid is found when both induced and chiral NNN forces are included 

(NN+NNN curve) . It should be noticed that these results are still preliminary, as 
not all relevant excitations of the ^He nucleus have been taken into account yet. 
In particular, the resonances are sensitive to the inclusion of the first six excited 
states of the ''He [2] (here only the g.s. is included). We will complete this study 
in the coming months. 

We have also obtained first results for the d-|-^He scattering phase shifts and 
the ground state of ®Li including both SRG-induced and chiral NNN forces. The 
model spaces adopted so far contain only the g.s. of the "^He and d nuclei within a 
A^max = 8 HO basis size. While calculations for larger A^max values and including 
pseudo-excited states of the deuteron (fundamental to model the deformation and 
virtual breakup of this nucleus) are under way, these preliminary results are already 
very promising. In particular, the inclusion of the NNN force produces a change 
in position and splitting of the '^Di and scattering phase shifts, which were 
not well described in our former calculation with only the A'A'' part of the SRG 
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interaction. The NNN force has also the effect of increasing the binding energy of 
the ®Li nucleus. 



4.2 Three-cluster dynamics 

A proper description of Borromean halo nuclei and three-body breakup reactions 
(but also virtual breakup effects) within the NCSM/RGM approach requires the 
inclusion of three-cluster channel states and the treatment of the three-body dy- 
namics. At present we have completed the development of the formalism for the 

Table 1 Ground-state energies of the ''■'^He nuclei in MeV. Both NCSM/RGM and NCSM 
calculations were performed with the SRG-N'^LO NN potential with A = 1.5 fm~^, and 
hf2 = 16 MeV HO frequency. Extrapolations were performed with an exponential fit. 
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treatment of three-cluster systems formed by two separate nucleons in relative 
motion with respect to a nucleus of mass number A ~ 2. Preliminary results for 
the g.s. energy of ^He within a ''He(g.s.)-|-n -|- n cluster basis and an A^max = 12, 
Tif2 = 16 MeV HO model space, are compared to NCSM calculations in Table [1] 
The interaction adopted is the SRG-N^LO A^A^ potential with yl = 1.5 fm"^ With 
such a low value of A, at A^max ~ 12 the binding energy calculations are close to 
convergence in both NCSM/RGM and NCSM approaches. The ~ 1 MeV difference 
observed is due to the excitations of the ^He core, included only in the NCSM at 
present. Contrary to the NCSM, in the NCSM/RGM the ^He(g.s.)+n -|- n wave 
functions present the appropriate asymptotic behavior. This will be essential in 
describing ^He excited states in the continuum, such as, e.g. the 1~ soft dipole 
resonance. Work towards the solution of the three-cluster NCSM/RGM equations 
with continuum boundary conditions is under way. 

5 Conclusions 

We presented an outline of the NCSM/RGM, an ab initio many-body approach 
capable of describing simultaneously both bound and scattering states in light 
nuclei, by complementing the RGM with the use of realistic interactions, and a 
microscopic and consistent description of the nucleon clusters, obtained via the ab 
initio NCSM. We discussed applications to fusion reactions that power stars and 
Earth based fusion facilities, such as the ^Be(p, 7)^6 radiative capture, and the 
■^11(0!, n)^He and '^He(ci,p)^He fusion reactions. Finally, we outlined the progress 
toward the inclusion of the NNN force into the formalism and the treatment of 
three-cluster dynamics, and presented an initial assessment of N N N-ioice effects 
on ^He(n,n)^He scattering, as well as preliminary results for the g.s. energy of 
^He computed within a '^He(g.s.)-|-n -|- n NCSM/RGM three-cluster basis. Since 
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the publication of the first results [ll[2l[5], obtained for nucleon-nucleus collisions, 
the NSCM/RGM has grown into a powerful approach for the description of binary 
reactions starting from realistic NN forces. A truly accurate ab initio description 
of light-ion fusion reactions and light exotic nuclei that encompasses the full NNN 
force and the three-cluster dynamics is now within reach. 
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